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ABSTRACT 

Cygnus X-2 is one of the brightest and longest known X-ray sources. We 
present high resolution optical spectroscopy of CygX-2 obtained over 4 years 
which gives an improved mass function of 0.69 ± 0.03 M Q (la). In addition, we 
resolve the rotationally broadened absorption features of the secondary star for 
the first time, deriving a rotation speed of i>sini = 34.2 ± 2.5 kms -1 (la) which 
leads to a mass ratio of q — M c /M x = 0.34 ± 0.04 (la, assuming a tidally-locked 
and Roche lobe-filling secondary). Hence with the lack of X-ray eclipses (i.e. 
% ^ 73°) we can set firm 95% confidence lower limits to the neutron star mass 
of M x > 1.27 M© and to the companion star mass of M c > 0.39 M . However, 
by additionally requiring that the companion must exceed 0.75 M Q (as required 
theoretically to produce a steady low-mass X-ray binary), then M x > 1.88 M 
and i < 61° (95% confidence lower and upper limit, respectively), thereby 
making CygX-2 the highest mass neutron star measured to date. If confirmed 
this would set significant constraints on the equation of state of nuclear matter. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual 
(Cygnus X-2, V1341 Cyg) — X-rays: stars 
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1. Introduction 

The distribution of neutron star masses provide fundamental constraints on the 
equation of state of condensed matter. Very precise determinations (e.g., Thorsett et al. 



1993| , Nice, Sayer, & Taylor |1996|) are available from time delays in millisecond radio-pulsars 
and are all consistent with 1.38 ± 0.07 M Q . Dynamical masses can also be obtained 
for accreting neutron stars in X-ray binaries. In particular, high mass X-ray binaries 
(HMXBs) would appear to be ideal candidates because the orbit of the two components 
are measurable through spectroscopic doppler shifts and X-ray pulses. This, combined 
with a determination of the inclination angle through X-ray eclipses (when favourable) 
leads to a full solution for the system parameters. Following this strategy, six neutron star 
masses have been obtained, and they all lie in the range 1.0-1.9 M (see van Kerkwijk, van 
Paradijs & Zuiderwijk |1995| and references therein). However, the uncertainties involved in 
these determinations can be quite large due to non-Keplerian perturbations in the radial 
velocity curves. These are very difficult to assess and are caused by a variety of effects such 
as stellar wind contamination, tidal distortion of the companion and X-ray heating. On the 
other hand, mass determinations in low-mass X-ray binaries (LMXBs) are very difficult to 
obtain both because their neutron stars do not (usually) pulse and the optical companions 
are normally overwhelmed by X-ray reprocessed radiation (van Paradijs & McClintock 



1995|) . Only in a few exceptional cases can the companion be detected (when it is evolved 



or during X-ray off-states) and thus it becomes feasible to extract dynamical information 
and set constraints on the system parameters. 

Cygnus X-2 is one of the few LMXBs in which the spectrum of the non-degenerate star 
is visible, contributing about 50 percent of the total visual flux. Estimates of the absolute 



magnitude of the donor and analysis of interstellar reddening (McClintock et al. 1984 ) 



imply a distance of ~8kpc and hence an X-ray luminosity of L x ~ 10 38 ergs . This large 
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luminosity is consistent with near-Eddington accretion rates onto a neutron star as typically 
observed in LMXBs (e.g., Hasinger et al. |1990|) ; the neutron star's presence in CygX-2 



is also indicated by the observation of X-ray bursts (e.g., Kuulkers, van der Klis, & van 
Paradijs |1995|) . The X-ray intensity and energy distribution are highly variable on different 
time scales, tracing out a £ Z' shaped track in the so-called X-ray colour-colour diagram 
with three distinct spectral states (Kuulkers, van der Klis, & Vaughan |1996| , and references 



therein). Hence it is classified as a Z source (Hasinger & van der Klis |1989| ), the variations 



of which are believed to be triggered by mass transfer rate changes. Multi-wavelength 
observations indicate that the strength of the UV continuum and the high-excitation lines 
are correlated with the states of the 'Z' diagram (Vrtilek et al. |1990| ; van Paradijs et al. 
1990| ). Spectral type variations (in the range A5-F2) with orbital phase have been reported, 



with the earliest spectral type occurring when viewing the X-ray irradiated hemisphere of 



the companion (Cowley, Crampton, & Hutchings |1979| ). 



From 1993 onwards we have collected high-resolution spectroscopy of Cyg X-2 with the 
aim of improving the system parameters and to resolve the rotation speed of the companion, 
thereby significantly refining the mass determination of the two components. In addition, 
we have searched for the presence of Li in the atmosphere of the companion star (see 
e.g., Martin et al. |1994j), the results and implications of which will be presented elsewhere 



(Martin et al., in preparation). 



2. Observations 

We obtained 40 red spectra (AA6340-6800, 0.40 A pixel" 1 dispersion) of Cyg X-2 at 
the Observatorio del Roque de los Muchachos using the 4.2m William Herschel Telescope 
(WHT), equipped with the ISIS triple spectrograph (Clegg et al. |1992|) , on the nights of 



1993 Dec 16-19, 1994 Oct 23-24, 1994 Dec 25, 1996 Aug 5, 1996 Dec 3, and 1997 Aug 1-7. 
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A 0.8-1.3 arcsec slit was used, depending on seeing conditions, giving spectral resolutions 
of 25-37 km s _1 . Cu-Ne arc spectra were obtained after every 1800 s exposure of the target. 
For the sake of the spectral classification and rotational broadening analyses we also 
observed a grid of 34 template stars using exactly the same spectral configuration (with the 
narrowest 0.8 arcsec slit) as for CygX-2. These stars cover a range of spectral types from 
AO to F8 in luminosity classes III, IV and V. 



3. Results 

Individual radial velocities were extracted through cross-correlation of the red spectra 
with the template star HR2489 (A9III), after masking out the broad Ha and He I A6678 
emission lines. A subsequent sine wave fit to the velocity points (Fig. 1) gave the following 
parameters (after renormalising the minimum reduced xl to 1): P = 9.8444 ± 0.0003 d; 
7 = -209.6 ± 0.8 kms" 1 ; K c = 88.0 ± likms" 1 ; T (HJD) = 2449339.50 ± 0.03, where 
To corresponds to the standard zero phase definition, i.e., inferior conjunction of the 
secondary star. These and all subsequent errors quoted are ±l-cr. In order to explore any 
non-symmetric effects in the radial velocity curve, e.g., artificial eccentricity induced by 



heating of the inner face of the companion (e.g., Davey & Smith |1992|) , we allowed the 
eccentricity e to be a free parameter in the fit. The presence of the eccentricity is only 
significant at the ~75% confidence level; in this case we get e = 0.024 ± 0.015. Therefore, 
we conclude that a circular orbit represents the best description of the data points and thus 
we assume that our measured K c corresponds to the true velocity semi-amplitude of the 
companion star. Our parameters are entirely consistent and substantially more accurate 
than those derived by Cowley et al. ( |1979| ) and Crampton & Cowley ( |1980| ). 

Combining K c and P in the expression for the mass function gives 



where q = M c /M x is the system mass ratio. 

Using the above ephemeris we obtained a doppler-corrected spectrum of CygX-2, in 
the rest frame of the secondary. A spectral type classification of A9 ± 2 for the companion 
was then derived through two different techniques: optimal subtraction of spectral type 
standards in the regions AA6380-6520, AA6620-6665, AA6700-6760 (further details to be 
found in Casares et al. |1997| ) and direct comparison of the Fel line ratio A6463/A6457. The 
former method was also applied to Doppler sums at the two conjunction phases (—0.05 to 
0.05 and 0.45 to 0.55), but no spectral type variation could be found (Fig. 2). Our result is 
in contradiction with Cowley et al. (1979), who claim orbital variations of the companion's 
spectral type (due to X-ray heating) in the range A5-F2. However, we note that their result 
is based on the ratio Call K to H, which is not a good diagnostic because the Balmer 
series (and perhaps also Call K) are clearly filled in by variable emission cores. Indeed, it 
has been noted (Kristian, Sandage, & Westphal |1967t Cowley et al. |1979|) that the Balmer 
lines in CygX-2 appear abnormally broad, which is expected if the absorption cores are 
filled in with emission. Therefore, we give more weight to the metal line ratios, which do 
not support spectral variations larger than two subtypes between the two conjunctions. We 
note that we actually detect an enhancement of the He I A6678 absorption at phase 0.5, 
together with an overall weakening of the metallic absorptions, produced by a ~60 percent 
increase in the continuum. These are indications of heating effects, although spectral type 
variations are not significant since the relative depth of the metallic lines is maintained 
throughout the orbit. We also note that the absence of spectral type variations is not due 
to changes in the overall average X-ray luminosity over the last ~20yrs. In fact, recent 
Rossi X-Ray Timing Explorer (RXTE) All Sky Monitor measurements (Wijnands, Kuulkers 
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& Smale |1996|) indicate the overall X-ray luminosity to be comparable to that measured 
during the observations of Cowley et al. (1979). 

In order to measure the rotational velocity of the companion star (v sini) only the 
highest resolution spectra (with a 0.8 arcsec slit) were employed. These correspond to 
the nights of 1996 Dec 3 and 1997 Aug 1-7. The technique consists of performing a x 2 
test on the residuals after subtracting different broadened versions of our templates from 
the doppler corrected sum of CygX-2. The template spectra were broadened through 
convolution with the rotational profile of Gray ( |1976|) which assumes a linearized limb 
darkening coefficient (e = 0.5 at 6500A and T e g = 7000 K, see Wade & Rucinski |1985| ). The 
broadened templates were also multiplied by a variable factor to account for the continuum 
excess of the accretion disc, prior to subtraction. Further details of this procedure can be 
found in, e.g., Marsh, Robinson, & Wood ( |1994| ). Most of our templates are significantly 
broader than CygX-2 and therefore only a subset of 20 templates (those with intrinsic 
broadenings ^15 kms -1 ) were considered in this analysis. These all give v sini values for 
Cyg X-2 in the range 30-39 kms -1 with a mean of ~35.4kms _1 . The use of the most 
appropriate template (A9 III) provides vsini = 34.2 ± 1.5 kms -1 (see Fig. 3.) The l-a 
uncertainty in vsini was derived by forcing the minimum \ 2 to increase by one (Lampton, 
Margon & Bowyer |1976j ), after renormalising the minimum reduced xt to 1. However, we 
note that this error is purely formal and does not include sources of systematic error such 
as the uncertainty in the limb-darkening coefficient. More realistic estimates would be given 
by varying e between its extreme values of and 1. By doing this we find a more realistic 
vsini = 34.2 ± 2.5 kms -1 which will be the value adopted by us hereafter. Because the 
donor star fills its Roche lobe and is synchronized with the binary motion, the rotational 
broadening provides a direct measurement of the binary mass ratio q through the expression 
(e.g., Home, Wade & Szkody [19861) : 
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49g 2/3 

,sinz = K c (1 + q) (2) 

Substituting our values of v sinz and iT c in Eq. 2 we find g = 0.34 ± 0.04 which, 
combined with Eq. 1 gives M x sin 3 z = 1.25 ± 0.09 M . On the other hand, the absence of 
X-ray eclipses provides a severe upper limit to the inclination of i ^ 73°. This provides 
the following 95% confidence lower limits to the masses of the components in CygX-2: 
M x > 1.27M and M c > 0.39 M . The results are plotted in Fig. 4. 



4. Discussion 

Our conservative lower limit on M c is well below the minimum secondary mass of 
~O.75M required in recent theoretical predictions for steady LMXB sources (King et al. 
1997] ) . Imposing this condition of M c > 0.75 M now yields a 95% confidence lower limit 



of M x > 1.88 M and a 95% confidence upper limit of % < 61°. This would make CygX-2 
the heaviest neutron star mass measured to date, and thereby provide support for stiff 
equations of state for nuclear matter (e.g., Cook, Shapiro, & Teukolsky |1994|) and would 



contradict the "softer" equations of state as described in, e.g., Brown & Bethe (|1994|) . On 



the other hand, if we assume a maximum possible mass of the neutron star of ~3.2M 
(e.g., Rhoades & Ruffini |1974| ), we infer a 95% confidence upper limit of M c < 1.28 M and 



a 95% confidence lower limit of i > 45°. We note that the inclination constraints are more 
or less consistent with recent ellipsoidal model fits to a compilation of BV photometric light 
curves (Orosz & Kuulkers, in preparation). 

The RXTE has discovered a maximum kHz QPO frequency at 1066-1171 Hz in 8 
persistent LMXBs which, if interpreted as the orbital frequency of the last marginally 
stable orbit, implies neutron star masses of 2.0 ± 0.2 M (Zhang, Strohmayer, & Swank 



- 9- 



1997 , see also Kaaret, Ford, & Chen 1997|) ). Recently kHz QPO have also been discovered 



in CygX-2 (Wijnands et al. |1997| ), so our mass estimate is in excellent agreement with 
that expected by Zhang et al. ( |1997| ). As they noted, this would be consistent with current 



evolutionary scenarios for LMXBs, where neutron stars would be born at 1.4 M but accrete 
at near-Eddington rates for ~10 8 years (van den Heuvel & Bitzaraki |1995| ). Dynamical mass 



determinations of other persistent LMXBs with evolved secondaries (e.g., GX1+4) will help 
to construct the distribution of neutron star masses and thereby allow new constraints to 
be set on the equation of state of nuclear-density matter. 



We thank Tom Marsh for the use of his optimal extraction routines and the MOLLY 
analysis package, and Tariq Shahbaz for discussions and assistance in preparing the figures. 
We also thank Vik Dhillon, Rene Rutten and Miriam Centurion for supporting the Service 
observations, and Ian Browne, Neal Jackson and Peter Wilkinson for an exchange of 
observing time. JC acknowledges support by the EU grant ERBFMBI CT961756. The 
WHT is operated on the island of La Palma by the Royal Greenwich Observatory in 
the Spanish Observatorio del Roque de Los Muchachos of the Instituto de Astrofisica de 
Canarias. 



-10- 

REFERENCES 

Allen, C.W. 1973, in Astrophysical Quantities, Athlone Press 
Brown G.E., & Bethe, H.A. 1994, ApJ, 423, 659 

Casares, J., Martin, E.L., Charles, P.A., Molaro, P. & Rebolo, R. 1997, NewA, 1(4), 299 

Clegg, R.E.S., Carter, D., Charles, P.A., Dick, J.S.B., Jenkins, C.R., King, D.L., & Laing, 
R.A. 1992, ISIS Astronomers' Guide, ING User Manual No. XXII 

Cook, G.B., Shapiro, S.L., & Teukolsky, S.A. 1994, ApJ, 424, 823 

Cowley, A.P, Crampton, D., & Hutchings, J.B. 1979, ApJ, 231, 539 

Crampton, D., & Cowley, A.P. 1980, PASP, 92, 147 

Davey, S., & Smith, R.C. 1992, MNRAS, 257, 476 

Giacconi, R., et al. 1967, ApJ, 148, L129 

Gray, D.F., 1976, in The Observations and Analysis of Stellar Photospheres. Wiley- 
Interscience, New York, p. 398 

Hasinger, G., & van der Klis, M. 1989, A&A, 225, 7 

Hasinger, G., van der Klis, M., Ebisawa, K., Dotani, T., & Mitsuda, K. 1990, A&A, 235, 
131 

Home, K., Wade, R.A., & Szkody P. 1986, MNRAS, 219, 791 
Kaaret, P., Ford, E.C., & Chen, K. 1997, ApJ, 480, L27 
King, A.R., Frank, J., Kolb, U., & Ritter, H. 1997, ApJ, 484, 844 
Kristian, J., Sandage, A., & Westphal, J.A. 1967, ApJ, 150, L99 



Kuulkers, E., van der Klis, M., & van Paradijs, J. 1995, ApJ, 450, 748 

Kuulkers, E., van der Klis, M., & Vaughan, B.A. 1996, A&A, 311, 197 

Lampton, M., Margon, B., & Bowyer, S. 1976, ApJ, 208, 177 

Marsh, T.R., Robinson, E.L., & Wood, J.H. 1994, MNRAS, 266, 137 

Martin, E.L., Rebolo, R., Casares, J., & Charles, P.A. 1994, ApJ, 435, 791 

McClintock, J.E., Remillard, R.A., Petro, L.D., Hammerschlag-Hensberge, G., & Proffitt, 
C.R. 1984, ApJ, 283, 794 

Nice, D.J., Sayer, R.W., & Taylor, J.H. 1996, ApJ, 466, L87 

Rhoades, C.E., & Ruffini, R. 1974, Phys. Rev. Lett., 32, 324 

Thorsett, S.E., Arzoumanian, Z., McKinnon, M.M., Taylor, J.H. 1993, ApJ, 405, L29 

van den Heuvel, E.P.J. , & Bitzaraki, O. 1995, A&A, 297, L41 

van Kerkwijk, M.H., van Paradijs, J., & Zuiderwijk, E.J. 1995, A&A, 303, 497 

van Paradijs, J., & McClintock, J.E. 1995, in X-Ray Binaries, eds., Lewin, W. H.G., van 
Paradijs, J., & van den Heuvel E.P.J. , Cambridge University Press, p. 58 

van Paradijs, J., Allington-Smith, J., Callanan, P., Charles, P.A., Hassall, B.J.M., Machin, 
C, Mason, K.O., Naylor, T., & Smale, A.P 1990, A&A, 235, 156 

Vrtilek, S.D., Raymond, J.C., Garcia, M., Verbunt, F., Hasinger, G., & Kurster, M. 1990, 
A&A, 235, 162 

Wade, R.A., Rucinski, S.M. 1985, A&AS, 60, 471 



- 12 - 



Wijnands, R., Homan, J., van der Klis, M., Kuulkers, E., van Paradijs, J., Lewin, W.H.G., 
Lamb, F.K., Psaltis, D., & Vaughan, B., ApJL, in press ( |astro-ph / 971 1219|) 

Wijnands, R., Kuulkers, E., & Smale, A.P. 1996, ApJ, 473, L45 

Zhang, W., Strohmayer, T.E., & Swank, J.H. 1997, ApJ, 482, L167 



This manuscript was prepared with the AAS IAT^X macros v4.0. 



-13- 




Fig. 1. — Radial velocity curve of the secondary star in CygX-2 folded on the ephemeris 
given in the text, plus the corresponding best sinusoidal fit. 
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Fig. 2. — Reduced y 2 of the optimal subtraction technique as a function of the spectral type 
of the templates. Different symbols indicate different luminosity classes: filled triangles for 
giants, open squares for subgiants and filled circles for main dwarfs. 
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Fig. 3. — Normalised spectra of CygX-2 and the A9III template HR2489 at the two 
conjunction phases. From top to bottom we plot the residual CygX-2 spectrum after 
subtraction of the broadened template, the doppler corrected sum of CygX-2 and the 
HR2489 spectrum broadened by 34kms _1 . 
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Fig. 4. — The mass of the companion (M c ) versus the mass of the compact object (M x ). 
Values for M c and M x for CygX-2 lie between the two solid lines which indicate the \-o 
error in the mass ratio q. The curved dashed line represents the constraint on the binary 
inclination by the lack of X-ray eclipes (i ^ 73°). The vertical dashed line indicates the 
maximum possible mass for the neutron star (3.2 M , e.g., Rhoades & Ruffini 1974). The 
horizontal dashed line indicates the mass of a main sequence star (1.7M , Allen 1973) of 
the same spectral type as observed for the companion in CygX-2 (A9). 



